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Abstract - We consider percolation on interdependent locally treelike networks, recently intro- 
duced by Buldyrev et al., Nature 464, 1025 (2010), and demonstrate that the problem can be 
simplified conceptually by deleting all references to cascades of failures. Such cascades do exist, 
but their explicit treatment just complicates the theory - which is a straightforward extension of 
the usual epidemic spreading theory on a single network. Our method has the added benefits that 
it is directly formulated in terms of an order parameter and its modular structure can be easily 
extended to other problems, e.g. to any number of interdependent networks, or to networks with 
dependency links. 
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On September 28, 2003, Italy experienced its most se- 
vere blackout in over 20 years pQ. The blackout's sever- 
ity was later attributed to the fact that an initial fail- 
le ure on the physical power grid disrupted not only the 

T^j- grid, itself, but also a computer network that depended 
■^j - on this grid for electricity [2(2] ■ Since the grid's substa- 
\l tions were, in turn, dependent on this computer network 
Q>for their regulation, further failures in the grid ensued as 
<^> communication among the stations was lost. Ultimately, 
,— ' recursive cascading failures throughout both networks oc- 
' . I curred, and both networks changed from percolating to 
non-percolating [IHI]- While this is a spectacular exam- 

■ pie of percolation on interdependent networks, it is, by no 

■ means the only one: in fact, most real- world networks can 
C3be seen as having some interdependency [5HZ]. A clear 

description of how perturbations propagate through such 
networks - i.e. of how perturbations can effect percolation 
(fragmentation) - is essential to understanding systems in- 
volving interdependence, including economic markets, in- 
terrelated technological and infrastructural systems, social 
networks, disease dynamics, or human physiology. 

Recent models of percolation on interdependent net- 
works have been described in terms of failures cascading 
back and forth between the networks [2HI]. While there 
is no doubt that percolation on interdependent networks 
can be seen as a cascading phenomenon, the mathemat- 



- r— I 

X 



ics behind such a description is cumbersome and far from 
transparent. Here we simplify matters by omitting all as- 
pects of cascading and treat percolation on interdepen- 
dent networks as an epidemic spreading process in com- 
plete analogy to ordinary percolation. If one wants to con- 
sider cascades explicitly, this can be done in a second step, 
after the phase transition itself is well understood. We 
also show that percolation on dependency networks [H|5] , 
which are single networks composed of both connectivity 
links and dependency links, can be described within the 
same paradigm. 

Both the theory of 2 4 and the present paper deal only 
with locally treelike random networks, for which mean field 
theory based on generating functions becomes exact in 
the large system limit. The cascade-based studies in PH2] 
considered site percolation networks from which a certain 
fraction 1 — p of nodes had been removed. On general 
networks (including lattices), this site dilution can lead 
to a topological modification of the networks j5J, but in 
the present cases it just leads to a trivial rescaling of the 
number of nodes and links. Thus we can restrict our- 
selves, without loss of generality, to the case p = 1 and 
only consider bond percolation, which further simplifies 
the discussion. 

In the following, we first recall the theory of epidemic 
spreading (ordinary percolation theory) on single networks 
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[TUHH], and then demonstrate how this can be easily 
adjusted to accommodate percolation on interdependent 
networks or on dependency networks. We only consider 
the limit of large networks, where the number of nodes 
N -> oo. 

(i) Single network - First consider a single (isolated) 
random network with mean degree z, with an epidemic 
spreading from some starting node. Let Si be the prob- 
ability that node i is infected during this epidemic (i.e. 
that it is part of the infinite percolating cluster). Its av- 
erage over all nodes, denoted as S, is taken as the order 
parameter of the model. The probability that node i is not 
infected is equal to the chance that none of its neighbors 
are infected through their remaining links: 

l-s^na-^). (!) 

m 

where the product runs over all neighbors of i and S'j is 
the probability that node j is infected through a randomly 
chosen edge not attached to i. When the graph is locally 
treelike, all S'j are independent. Averaging Eq. ((T|) over 
all nodes gives then 

5=1- - S') k = 1 - G (l - 5'), (2) 

k 

where p(k) is the probability that a node has k links and 
Go(x) is the corresponding generating function Gq(x) = 
^2 k p(k)x k . Similarly, one can write down the equation 
for 5' 

S' = 1-Y = i_ Gl (i_5') ( 3 ) 

k 

where z = G' (l) and G\(x) = G' (x)/z. For any given 
degree distribution we can first solve Eq. (|3) to obtain S', 
and then insert it into Eq. (J2J to obtain S. The percolation 
transition threshold z c is given by S = for z < z c and 

s > o for z > Zc main]. 

For instance, the degree distribution of Erdos-Rcnyi 
(ER) graphs [HJH2] * s Poissonian. Therefore, 

G (l-S') = J2^^-S') k = e- zS ' 

k 

= Gx(l-S'). (4) 

Thus, S = S' and Eqs. © and © give simply 5 = 1- 
e~ zS . Defining /(5) = 5 — 1 + e~ zS , one can find the 
solution S(z) by solving f(S) = graphically, as shown 
in Fig. [TJ A continuous ('second order') phase transition 
is clearly evident in the inset of Fig. [T] 

Equations ([2]) and (|3|), are the order parameters for a 
single network, where 5 represents the probability that 
a randomly chosen node places in the infinite percolating 
cluster and 5' is the same probability, but defined when we 
pick an edge randomly and look at the end node. These 




Fig. 1: Graphical solutions for ordinary percolation on a single 
ER network. The inset shows the continuous change of order 
parameter S as z increases. 

equations act as fundamental 'building blocks' or 'mod- 
ules' for treating analogously the probability, on fully or 
partially interdependent networks, to be connected to the 
infinite percolating clusters. 

(ii) Two fully interdependent networks - Consider now 
two networks A and B, where each node in A depends 
only on one node in B and vice versa. In order for a 
node in network A to be part of the percolating cluster, 
its partner in B must also be part of that cluster. Since 
this mapping is one-to-one we can merge each node in 
A with its partner in B to have one set of nodes, each 
with two sets of links. We define .A£>-clusters as subsets 
of nodes connected both in A and in B. More precisely a 
set of nodes C = ■ ■ • 4} is an _4£>-cluster if any two 

points i,j E C are connected by two paths: one path using 
only links 6 A and nodes only £ C, and the other using 
only links e B and also using nodes only £ C. We do not 
allow paths that involve nodes outside C, so ^IB-clusters 
are self-sustaining 

The probability that any node belongs to the infinite 
.4£>-cluster is equal to the probability to be linked to it 
both via A- and via S-links. That is, a node looks out at 
its ^4-links to see if it has a neighbor on the percolating 
.4£>-cluster. It also looks out via its B links. Only if it 
has a neighbor via both sets of links is it a member of this 
cluster. Therefore 5 is simply a product of the right hand 
side of Eq. © for networks A and B, 

S=(1-G£(1-S' A ))(1-G$(1-S' B )), (5) 

where the superscripts (subscripts) A and B refer to net- 
works A and B respectively and Sa = 5g = 5. Here, S' A 
(and analogously S' B ) is defined as the probability that a 
node reached by following a random ^l-link is in the AB- 
cluster. For this to happen, its partner node - which is a 
random node from the point of view of network B - has 
also to be connected to the „4£>-clustcr via £>-links. S' A 
and S' B arc different from each other since they depend 
on the degree distribution of each network. When choos- 
ing edges at random, the end node of a randomly chosen 
edge in network A belongs to the ^IB-cluster only when 
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Fig. 2: Graphical solutions for the percolating cluster on two 
fully coupled ER networks with the same mean degree z. The 
inset shows a discontinuous change of order parameter S as z 
increases. 

its partner node in network B is concurrently a member 
of this cluster, and vice versa. The probabilities of these 
events occurring - S' A and S' B - are given by 

S' A = (l-G?(l-S' A ))(l-G$(l-S' B )), 

S' B = (1-G?(1-S' B ))(1-G$(1-S' A )). (6) 

Using this and Eq. @ for two interdependent ER net- 
works with mean degrees z A and zg gives 

S= (1 -e-^ s )(l -e- ZflS ). (7) 

In particular, when z A = z B = z, the order parameter 
obeys 

S= (1 -e~ zS f. (8) 

Defining g(S) = S — (1 — e~ zS ) 2 , one can find graphically 
the value of S as a function of z that solves g(S) = 
(see Fig. [2]). This solution shows a discontinuous ('first 
order') phase transition (inset of Fig. [2]), contrary to the 
result for the single network in Fig. [TJ Demanding g(S) = 
g'(S) = we find a critical point S c = 0.511699 • • • and 
z c = 2.455407 • • • . The value of z c agrees with that given 
for ER networks in [2 [9] and our explicitly determined 
value of the order parameter just above the critical point, 
S c , can also be obtained from appropriate combinations 
of their results. 

For the general case ^ zg, defining 

h(S) = S- (1 -e- ZAS )(l -e" 2BS ), (9) 

we obtain a line of discontinuous transition points from the 
conditions h(S) = and h'(S) = 0. Assuming h(S) = 0, 
the second condition can be written as 

h'(S) = 1 - z A e- ZAS (l - e- ZBS ) - z B e~ ZBS {l - e~ ZAS ) 
z A S e- ZAS z B S e- zeS 
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Fig. 3: (a) S as a function of z A and zb for Eq. (0). (b) The 
discontinuous phase transition line on the plane of z A and zb- 
It is symmetric about the line z A — zb- In the inset, the jump 
size S c at the transition line is shown as a function of the ratio 
za/zb- 

where x = z A S and y = z B S. This gives a one-parameter 
set of solutions y(x), from which S c can be obtained using 
Eq. (J7J). Finally, z A and zb are obtained by z A = x/S and 
Z B = y/S. The resulting discontinuous phase transition 
line is shown in Fig. [3]^b) , where z A = 1 and zg = 1 act 
as asymptotes. In the inset of Fig. (3Jb), the jump size S c 
is shown as a function of the ratio z A j 1 zb (only z A j zb < 1 
is shown, since the transition line is symmetric about the 
line z A = zb). The full dependence of S on z A and zb is 
shown in Fig. EHa). 

(iii) An arbitrary number of interdependent networks - 
Our approach is easily extended to treat coupling of more 
than two networks [13] ■ For M networks, Eqs. ([5J) and ^ 
can be simply replaced by 

M 

m— 1 

c* _ — Gy(i — s f m ) , , 

This is due to the fact that in order to be part of the infi- 
nite percolating cluster each node must, by definition, be 
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connected to it via all of the M networks. Order parame- 
ters and transition points can be obtained in an analogous 
way and the transition is discontinuous for any M > 1. 

(iv) Partially interdependent networks - Assume now 
that two networks are not totally interdependent. In net- 
work A only a fraction qA of the nodes are dependent on 
a node in network B. To be on the percolating cluster, a 
node has to be connected to the cluster via A links and 
the node on which it depends has to be connected to the 
cluster via B links. Similarly a fraction gg of nodes in net- 
work B depend on a node in A. Note that here one node 
from a network depends only on one node from the other 
network, i.e., each node can have only one dependency 
link, which can be unidirectional or bidirectional [3J. In 
that case, in general we must expect that the two order 
parameters S A and Sb arc different. They indicate the 
chance that a randomly picked node in A (resp. B) is a 
member of the percolating A6-cluster. 

Let us first discuss the symmetric case qA = qs = Q- If 
a given node in network A does not depend on network 
B (which happens with probability 1 — q), it is a mem- 
ber of the percolating ^IB-cluster iff at least one of its 
neighbors in A is also a member of that cluster. On the 
other hand, with probability q, the node does depend on 
a node in B. In that case, in order for it to be part of the 
percolating _4£>-cluster, its dependency partner also must 
be connected via network B to at least one node in that 
cluster. Therefore, the probability Sa can be expressed by 
summing the conditional probability to be connected mul- 
tiplied by the corresponding probability to be dependent 
(q) or not (1 — q); 



9(1-<#(1 

+(1-3)(1- 
(1 - Gjftl - 



-^))(1-G 6 (1 

■G$(1-S A )) 
S' A ))0—qG*(l 



S'b)) 



S'b))- (12) 



Similarly, 

S' A = (1 - Gf(l S' A ))(1 qG*(l S' B )). (13) 

By symmetry, another pair of conditions exists for network 
B: 

S B = (l-G*(l-S' B ))(l-qG£(l-S' A )), 

S' B = (l-G*(l-S' B ))(l-qG$(l-S' A )). (14) 



More generally, if q A qs, 



S A — 
Sb = 

Sr = 



(l-G£(l-S' A ))(l-q A G%(l-S' B )), 
(1 - G?(l - S' A ))(1 - QaG®(1 - S' B )), 
(l-G$(l-S' B ))(l-q B G$(l-S' A )), 
(l-Gf(l-S' B ))(l-q B G$(l-S' A )). 



(15) 



Again, these four coupled equations involve nothing more 
complicated than (weighted) products of the network- 
specific fundamental factors from Eqs. and ([3]). To- 
gether, they give the behavior of the order parameters S A 
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Fig. 4: (Color online) (a) Percolation transition points for two 
partially coupled ER networks with the same mean degree z at 
different values of the dependency q. (b) The jump size Sc at 
the transition point. The first order phase transition line (q > 
1/3) meets the second order phase transition line (q < 1/3) 
at q — 1/3. The red dotted line indicates the tricritical point 
(q a ,z a ) = (1/3,3/2). 



and Sb- These expressions arc completely equivalent to 
the more complicated results in [3J. If q A = q B = q = 1, 
these equations reduce to Eqs. and ([5]). 

If q A = q& = q ^ 1 and if the two networks have the 
same degree distribution, Sa = S b — S. For example, if 
two ER networks having the same mean degree z are cou- 
pled with the dependency probability g / 1, the solution 
is simply 

S = (1- e- zS )(l- qe' zS ). (16) 

The solution S can now be expressed as a function of both 
z and q. When 5 = 0, S(z) shows a continuous transition, 
since now both networks are fully independent. On the 
other hand, when q = 1, the transition is discontinuous. A 
crossover from second order behavior to first order behav- 
ior at the tricritical point q c is observed as q is increased 
from to 1. In analogy to Eq. ©, we now define 



h(S) = S - (1 - e- zb )(l - qe 



-zS\ 



(17) 



The tricritical point is found by demanding S = h(S) = 
h'(S) = h"{S) = 0, which gives q c = 1/3 and z c = 3/2 0. 

For q > q c , we can find the first order transition point 
by considering h(S) = and h'(S) = 0. The subsequent 
analysis is straightforward and follows closely our previous 
method. Its results are displayed in Fig. [4] 

Let us briefly discuss the case qA = Qb = 1 and za ^ z B 
of ER networks, which was not treated before. In this case, 
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Fig. 5: (Color online) Phase transition line at q = 1/2 for 
different 2u and zg. The black solid line indicates the first 
order transition line and the blue dash-dot line represents the 
second order transitions. They meet two tricritical points (red 
bullets). The tricritical points move along the red dotted line 
as q changes. When q = 1/3, the tricritical points meet at 
za = zb = 3/2. 

also S A ^ Sb , and we have to solve the coupled equations 

h A (S A ) =S A -(1- e-^)(l - qe- z ^) = 0, 
h B (S B ) = S B - (1 - er ZBS ' 3 ){l - qe- z * SA ) = 0. (18) 

The percolation transition is obtained by imposing in ad- 
dition h'^{S A ) = h' B (St3) = 0. Defining again x = z A S A 
and y = zgSg, we find 

f 1 - - _»_ > \ = x . (19) 

V e x - 1 J \ ev - I J e x - q ev - q 

For given q this is solved graphically (Fig. [5J . The values 
of S A and Sg at the transition point are then obtained 
from Eq. (fT5|) . If they vanish, the transition is continuous, 
otherwise it is discontinuous. Results are shown in Fig. [3] 
for q = 1/2. If q decreases, the two tricritical points move 
together. They coalesce at z A = z& = 3/2 when q — > 1/3. 

(v) Dependency networks - A dependency network [4j 
[5] is a single network that contains two types of links, 
connectivity links and dependency links. In the simplest 
case each node in the network depends on one other node 
in that network and all dependencies are mutual. In order 
for a node to be connected to the infinite self-sustaining 
cluster both it and its partner must be part of that cluster. 
For random networks this leads immediately to Eqs. ([5]) 
and © where the subscripts A and B are dropped since 
there is only one network. In the case that only a fraction 
q of nodes have dependency links we get Eqs. (fl"2f and 
(|13|) . again dropping the subscripts labelling the networks. 
Again the type of transition depends on the value of q, 
with a tricritical point separating the two regimes. 

In the most general case where a node has m depen- 
dency links with probability p(m) the resulting equations 



are 

JV-l 

S = ]Tp(m)(l-Go(l-S')r +1 , 

m=0 

JV-1 

S' = (1 - Gi(l - S')) £ p(m)(l - G (l - S')) m - 

Again the precise behavior reflects the extent to which the 
network is dependent, with a low dependency regime ex- 
hibiting a continuous transition in the universality class 
of ordinary percolation and a high dependency regime ex- 
hibiting a first order transition, with crossover controlled 
by a tricritical point. Similar arguments can be used to 
derive the general case for interdependent networks where 
a single node has dependencies to m other networks with 
probability p(m). In that case the equations are more 
complicated because each network's structure may be dif- 
ferent and they will have different order parameters, but 
the arguments used to derive the equations are precisely 
the same. The above does not apply to networks with 
directed (non-mutual) dependencies, for which different 
arguments apply. 

In summary, we consider percolation on various inter- 
dependent or dependency networks, pointing out the close 
analogy to epidemic spreading on single networks with- 
out these dependencies. Our arguments are much more 
straightforward, both conceptually and mathematically, 
than previous ones built on cascades of failures. We should 
however stress that they apply, like those of [2"j flH8"lfT3] only 
to random locally treelike graphs. For interdependent net- 
works that are correlated with each other Q3] or that are 
spatially embedded JSJUH] the transition is in general not 
first order, and the interdependency can make the tran- 
sition even less sharp than in ordinary percolation [91IT5]. 
For these more realistic cases no analytical theory is yet 
available. 
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